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The increasing transmission capacity needs in a future energy system raise the question how 
associated costs should be allocated to the users of a strengthened power grid. In contrast to 
straightforward oversimplified methods, a flow tracing based approach provides a fair and consistent 
nodal usage and thus cost assignment of transmission investments. This technique follows the power 
flow through the network and assigns the link capacity usage to the respective sources or sinks using 
a diffusion-like process, thus taking into account the underlying network structure and injection 
pattern. As a showcase, we apply power flow tracing to a simplified model of the European electricity 
grid with a high share of renewable wind and solar power generation, based on long-term weather 
and load data with an hourly temporal resolution. 
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I. INTRODUCTION 

The main drivers of the increasing integration of re¬ 
newable energy into present and future energy systems 
are climate change mitigation and sustainability T]. Due 
to its fluctuating and decentralized nature, the renew¬ 
able energies represent a challenge in terms of security 
of supply, robustness and cost efficiency. In this regard, 
detailed techno-economic models of the energy system 
are invaluable for the evaluation of well-defined scenarios, 
but often lack the possibility to assess fundamental mech¬ 
anisms and interrelations due to their sheer complexity 
and multitude of parameters. On the other hand, over¬ 
simplified abstract models often neglect fundamental cor¬ 
relations and structures present in real systems. At this 
point, the applied theoretical physics of complex renew¬ 
able energy networks represents an important bridge be¬ 
tween these poles I2H5I- The incorporation of suitably ag¬ 
gregated long-term weather-determined renewable gen¬ 
eration and historic load data accounts for the relevant 
spatio-temporal correlations, while the usage of simpli¬ 
fied network and dispatch models reduce the complex¬ 
ity of the problem and limits the number of parameters. 
Such an approach allows to identify fundamental inter¬ 
relations, develop and test new concepts - often inspired 
from physics - and thus is able to guide the research of 
more in-depth models. This point of view is shared by 
several recent contributions from the physics community 
addressing challenges for modern complex electricity in¬ 
frastructures, with a strong focus on power grid synchro¬ 
nization, stability, and control EMU- 

Following this line of research, in [TST - fTS! it has been 
shown that for a highly renewable European electricity 
network, the overall transmission capacity needs can be¬ 
come quite large. In this respect, a simple, but important 


question emerges: who is going to pay for this? Straight¬ 
forward propositions distribute the overall transmission 
capacity onto the nodes proportional to their average 
load, or count only the links l{n ) attached to a node n. 
Although the average load is a first approximation of the 
respective network usage, it does not take into account 
the geographical position of a node in the network, nor 
any correlations in the import/export distribution of the 
system. Power flow tracing incorporates such relation¬ 
ships and thus is a vital ingredient of a more appropriate 
usage measure underlying a “fair” and consistent cost 
allocation procedure. 

This technique, also called average participation, has 
been pioneered by m and HOI [2I]- It has been de¬ 
veloped and suggested for transmission congestion man¬ 
agement and pricing [22], for instance in the context of 
inter-Transmission System Operator (TSO) compensa¬ 
tions, that is the tariffication of cross-border flows be¬ 
tween different European grid zones [2,'ll . Applying the 
principle of proportional sharing, the tracing algorithm 
acts on a given power flow solution and provides a parti¬ 
tion of the total flow on the network into different partial 
flow patterns emerging from the respective sources, or al¬ 
ternatively ending in the respective sinks. These partial 
flow patterns can be represented by colours according 
to their originating (alternatively ending) node. Techni¬ 
cally, flow tracing then corresponds to colour diffusion on 
a directed flow graph embedded in the original network. 

In this paper we introduce a new link usage measure 
based on the flow tracing method and apply it to a sim¬ 
plified model of a highly renewable European electricity 
network. Using a notation based on propagating colour 
vectors, we reformulate the tracing technique using con¬ 
cepts familiar from diffusion processes on complex net¬ 
works m- Considering an ensemble of different flow 
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patterns on a given network, the flow tracing techniques 
yield an ensemble of coloured partial flow patterns. The 
new usage measure integrates this information into a sin¬ 
gle value for each source, sink and link in the network. 
The application of these methods to the European elec¬ 
tricity network then provides information about the re¬ 
spective grid usage of different countries, and the com¬ 
position of their power imports and exports. Using ap¬ 
propriate cost functions, this usage measure can be inte¬ 
grated into future cost allocation schemes which provide 
a “fair” pricing of network related costs to individual im¬ 
porters and exporters. 

This paper is organized in the following manner: Sec¬ 
tion [TT] provides a simplified description of renewable en¬ 
ergy networks. Flow tracing on a network is the topic of 


Section |III| Section IV applies flow tracing to a simpli¬ 
fied highly renewable European electricity network, and 
includes an ensemble averaging over fluctuating flow con¬ 
figurations. A conclusion and outlook is given in Sec¬ 
tion |v] 


II. RENEWABLE ELECTRICITY NETWORKS 


In a renewable electricity network, like the one illus¬ 
trated in Figure [l] each node n comes with a renewable 
power generation G^(t) and a load L n (t). Their temporal 
averages ( G%) = 7 n (i„) define the renewable penetra¬ 
tion 7 „. Even for 7 „ = 1 the renewable power generation 
at a specific time t almost certainly will not match the 
load. It requires a balancing infrastructure to deal with 
this mismatch: 

G*(t) - L n (t) = B n (t) + P n (t) . ( 1 ) 

The nodal balancing B n (t) describes the backup power 
generation G^(t) = — min{S n (t), 0} when B n (t) < 0, 
and the curtailment C n (t) = max{B n (t), 0} of excess 
power when B n (t) > 0. A positive injection P n (t) > 0 
allows the node to export power to other nodes in the 
network and to reduce its curtailment. In case of a neg¬ 
ative injection P n (t) < 0 the node is importing power 
from other nodes, which reduces its own backup power 
generation. The sum Pn(t) = 0 over all nodal in¬ 

jections has to be zero, representing flow conservation 
over the network with lossless power flows. N represents 
the number of nodes in the network. 

Equation ([I]) can be understood as an “actio = reactio” 
equation. The mismatch on the left-hand side drives the 
networked system. The right-hand side represents the 
response of the system. In principle, more response types 
can be added to the right-hand side, like storage and 
coupling to the heating and transportation sectors. 

There are various ways how to divide the mismatch on 
the left-hand side of equation 0 onto the balancing and 
the injection. A particularly simple and appealing inter¬ 
action between B n and P n is described by synchronized 


balancing [BJ: 


B n (t) 


E (Gft(t)--Mt)) 


\m—l 


(Ln) 

J2k(Lk) 


( 2 ) 


It distributes the total mismatch onto each node in pro¬ 
portion to its average load. Together with equation ([I]) 
the synchronized balancing ([ 2 ]) fixes the injection pat¬ 
tern P n (t ). The constraint J2n=i^(t) = 0 is fulfilled 
automatically. 

The injection pattern determines the active power 
flow Fi along the links l = {m—¥n) [25] : 


Pn(t) = YB nrn e m {t ), 

m 

Fm^n(t) = b mn (6 m (t) - 9 n (t)) . 


(3) 

(4) 


Here, the DC approximation has been evoked, where link 
resistances can be neglected |25l [26i. As in [6], for conve¬ 
nience all reactances are set to one and the dependence 
of the reactances on line length is disregarded, as this 
leads only to slightly modified results. The susceptance 
matrix B nm = (J2 k Kk)Km - Km then simply becomes 
identical to the Laplacian of the network graph, where 
Km = 1 if nodes n and m are directly connected by a 
link and b nm = 0 otherwise. The variables 9 n represent 
the voltage phase angles. 

Figure [l] illustrates the nodal injections and the result¬ 
ing flows on the links in a simplified European transmis¬ 
sion grid for one time instance. Data from a Renew¬ 
able Energy Atlas [2] have been used for the renewable 
power generation G^(t) = G^ ( t ) + Gjj(f), where spatio- 
temporal wind velocity and solar radiation fields have 
been converted into country-specific long-term wind and 
solar power generation time series with hourly resolution. 
Historical load data provides the demand or load L n (t) 
per country with the same temporal resolution [2]. For 
all the results shown here in this paper a fixed renewable 
penetration parameter 7 n = 1 and a fixed wind fraction 
<Gr>/(G«) = 0.8 have been used. 

Sampling over time produces backup power distribu¬ 
tions P n {G^) for the nodes and flow distributions Vi{F{) 
for the links. The averages (G^) of the nodal distri¬ 
butions describe the required average backup energy. 
High quantiles of these distributions describe the required 
backup capacities, which secure the power supply almost 
surely at all times. Backup energy and capacity have 
been discussed for example in [3 E], Here, we are in¬ 
terested in the required transmission link capacities JCf 
for unconstrained power flow through the network. They 
follow from the high quantiles of the distributions for the 
absolute link flows fi = |E]|: 

ricT(q) 

q= Wi)dfi. (5) 

Jo 

We adopt q = 0.99. Other definitions for the trans¬ 
mission link capacities have been given for example 
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FIG. 1: Injection pattern and resulting power flows on a sim¬ 
plified European transmission network nz). A typical winter 
night hour has been picked from the Renewable Energy Atlas 
developed in [3j. 


in nans]. The sum over all link capacities 

= ( 6 ) 

;=l 


also [2DH22]- Using the language of physics, flow trac¬ 
ing is related to a vector diffusion process on a directed 
flow graph embedded in a network. In the following we 
explain the details. 

Fi n (t) describes the part of the flow Pj(f) along link l, 
which has originated in the exporting node n. Their ratio 
is denoted as Ci n {t) = Fi n (t)/Fi{t) and forms one compo¬ 
nent of the colour vector cj(t) = (cn(t),... ,CiN(t)) T of 
link l ; here the superscript T stands for the transpose of 
the vector, and not for power transmission as it has been 
used for example in equation ([ 6 ]). The one-norm of this 
vector is equal to ||q||i = J2 n c in = !■ 

The colour vector along the directed link l = 
which has the same direction as the flow Fi ( t ), is identical 
to the colour vector of the upstream node to, cj(t) = 
Cm(t). With other words, the upstream node copies its 
colour vector onto all of its outgoing flows F m ^ n and 
(P m )_, where n are direct neighbours of to. Here we 
have used the notation ci for the colour vector on link Z, 
and c n for the colour vector on node n. Note that both 
types of vectors have N components denoted as Q m and 
c n m, respectively, which correspond to the part of the 
(out-)fiow originating in node to. 

The colour vector c m of node to. is determined by merg¬ 
ing all of its ingoing flows and their respective colour 
vectors: 


^ ' F k ^ m C k T (Pm) + 6m — F rn C rn , (7) 


determines the overall transmission capacity of the net¬ 
work. L represents the number of links in the network. 

Based on the data from the Renewable Energy At¬ 
las [ 2 ] the overall sum i> turns out to be eight times 
bigger than the overall current interconnector strenghts 
(Net Transfer Capacities) across Europe. For some of the 
links this ratio becomes even larger. Together with the 
average nodal loads ( L n ), the specific values for the trans¬ 
mission link capacities have already been listed in S3- 
Although a flow scheme different from ([2]) has been used 
there, the resulting transmission link capacities are qual¬ 
itatively similar to the present ones. The big increase 
in overall transmission capacities raises the question on 
who has to pay for this investment. A power flow tracing 
based network usage measure helps to clarify this point. 

III. FLOW TRACING 

For a grid usage measure, the link capacity K.J should 
be related to those nodes who are using the link Z for 
their own power exports and imports. The fraction JCf n 
associated with node n should be proportional to its link 
usage. Flow tracing allows to determine this link usage 
by either following where the nodal exports (P„(f))+ = 
max{P„(f), 0} are flowing to, or by tracing back where 
the nodal imports (P n (t))_ = — min{P„(f),0} are com¬ 
ing from. Power flow tracing has been first proposed 
in m using the principle of proportional sharing; see 


where 

*%= £ F k ^ m + (P m ) + ( 8 ) 

feeN™ 

represents the total in-flow into node to. The sum goes 
over all direct neighbours fc, which produce ingoing flows 
for to. e m is a unit vector with components e m i = 8 m i- 

This procedure represent a colour diffusion on a di¬ 
rected network without loops. Equation ([7]) mixes all 
colour vectors flowing into a node, and the mixed colour 
vector is then copied onto and distributed via all of the 
outgoing flows. 

The algorithmic solution of equation ([7]) for the nodal 
colour vectors uses a downstream procedure. It is ex¬ 
emplified in Figure (2jr+b. The first node n = 1 injects 
Pi = 2 units into the network. Since it has no ingo¬ 
ing flows on links, the only ingoing colour is the node’s 
own injection. This results in the colour vector ci = 
( 1 , 0 , 0 , 0 , 0 ) T , where the first component has been given 
the colour blue in Figure [!>■ Analogously, we get the 
colour vectors c *2 = ( 0 , 1 , 0 , 0 , 0 ) T and C 5 = ( 0 , 0 , 0 , 0 , 1 ) T 
for the second and fifth node, which are like the first node 
only pure upstream nodes with positive injections. The 
colours of the second and fifth component are illustrated 
in yellow and red, respectively. The colour vectors of 
the first, second and fifth node are copied onto the first, 
second and fourth link, respectively. The third node has 
ingoing flows along the first two links and no own positive 
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FIG. 2: Illustration of a simple five-node network with (a) 
an injection and a resulting flow pattern, (b) the flow tracing 
of exports, and (c) the flow tracing of imports. 


injection. This fixes the colour vector of the third node 
to c 3 = ( 1 / 2 , 1 / 2 , 0 , 0 , 0 ) T , which is shown as the nodal 
50% blue / 50% yellow mix in Figure §>• This colour 
vector is then copied onto link l = 3, which carries the 
link-based outflow of the third node. The node n = 4 
receives ingoing flows via the links l = 3 and 4, and has 
no positive nodal injection. The weighting of the two 
ingoing colour vectors with the respective flow strengths 
then yields the colour vector C 4 = (1/6,1/6,0,0, 2/3) T . 

Figure illustrates the export picture of flow trac¬ 
ing, where link usage is caused by exporting nodes. This 
can of course also be turned around to the import pic¬ 
ture of flow tracing. In this dual picture, the role of the 
sources and sinks in the injection pattern is exchanged, 
i.e. P n (t) —> —P n {t). It is illustrated in Figure [2J3, and 
allows to define link usage caused by importing nodes. 


IV. APPLICATION OF FLOW TRACING TO A 
SIMPLIFIED HIGHLY RENEWABLE 
EUROPEAN ELECTRICITY NETWORK 


The mismatch between the renewable power genera¬ 
tion and the load on the left-hand side of equation ([l]) 


changes with each time step. This causes also the nodal 
balancing B n (t), the nodal injection P„(t), the power 
flows Fi{t) on the links, and the colour vectors c n (t) and 
ci(t) of the nodes and links to depend on time. Figure [3] 
shows the time-sampled values for one colour compo¬ 
nent ci n , which have been obtained with the 8 years long 
hourly Renewable Energy Atlas data from 2j. About 
half of the time the values are ci n = 0 , which is because 
node n is not injecting (ejecting) in the export (import) 
picture. For the other half the values appear to scatter 
all over. The figur shows that in general there are various 
injection patterns P n (t ) which for a given sampling res¬ 
olution yield the same absolute flow fi(t) on a link, but 
different colour components ci n (t). The conditional aver¬ 
age (ci n \fi) over the respective “slices” corresponding to 
the same absolute link flow f) is also shown in Fig. [3] It 
shows a smooth dependence on /;, which is different for 
different n. Note however, that the nodal sum is always 
equal to ( c ln\fi) = 1- Based on this conditional 

average we will now derive an expression for a fair nodal 
assignment of link capacity usage. 

Imagine the link capacity Kj = f Q ‘ dK being built by 
increments dK,. Every increment dK = dK n is used 
by several nodes, and decomposes into the nodal usage 
contributions dK n = wi n (K)dK. The fraction Kf n of the 
overall link capacity Kf = Kj n used by node n then 
becomes 




wi n (K)dK . 


(9) 


As we go from capacity K to K + dK , all flows larger 
than K make use of this additional capacity incre¬ 
ment dK. Thus, the weight can be expressed as 


Wln{K)= P|(/l|/l>^)<Qn|/l)^J 

■Ik. 


where the conditional probability 


ViUilfi >K) = 


l-7f(/C) 


( 10 ) 


( 11 ) 


is proportional to the flow distribution Vi(fi), and 
Vf (K) = f Q Vi(fi)dfi represents the cumulative flow dis¬ 
tribution. Insertion of (10) and (11) into (|9j) leads to the 
expression 


r^i fjir r^i 

IC Tn = J o l ^ V C { K) J K WOfanl/lWl (12) 

for the fraction of link capacity used by node n. It 
is straightforward to check that this assignment fulfills 

Z n K? n = lCf. 

Based on the 8 years long hourly Renewable Energy 
Atlas data from [ 2 ], the resulting relative link capacities 
Kjn/Kj are summarized in Figure |4j Each row in this 
figure represents a link, and each column represents a 
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FIG. 3: Colour component Ci„(t) as a function of the power 
flow fi(t) = |Fl(f)|: (a) based on the flow tracing of exports, 
and (b) based on the flow tracing of imports. The link l = 
(CH,IT) between Switzerland and Italy, and the node n = DE 
(Germany) have been chosen. Each dot represents one hour 
contained in the 8 years long Renewable Energy Atlas data 
from [2|. The red curve represents the conditional average 
(cin\fi)- F° r the latter, the extreme flow events fi > Kj , 
which happen in 1 — q = 1% of the time, have been mapped 
onto /( = K.J. 


nodal country. The countries have been ordered accord¬ 
ing to their average load (L n ), with Germany having the 
largest load and Luxembourg the smallest. 

Germany is using all the links contained in the net¬ 
work. This is also illustrated in Figure UK The next 
biggest countries still use almost all of the links. The 
smaller a country becomes, the less number of links are 
used; see also Figure [HJd, which illustrates the case of 
small-medium sized Austria. Intuitively this is clear. 
Due to their big size, large countries are able to inject 
/ eject large exports / imports, which then penetrate the 


whole network. A quantitative measure for this penetra¬ 
tion is described by the average nodal export and import 
transfer functions: 


f 

< ~'n 

T 


i(t) = 
l (t) = 




< 0 Pn) + ) 

(^nm(Fn)-) 

((Pn)-) 


(13) 

(14) 


The former describes how much of the positive injection 
(P n )+ at node n is exported to a sink node to. The lat¬ 
ter expression describes how much of the negative injec¬ 
tion (P„)_ at node n is imported from a source node to. 
The export transfer function is also illustrated in Fig¬ 
ure [5] Whereas the German exports spread over the 
whole European network, smaller Austria mainly serves 
first neighbours and bigger second neighbours with its 
exports. The import transfer functions are not shown. 
They are similar to the export transfer functions. 

We now compare the flow tracing based usage mea¬ 
sure with two simplified propositions which already have 
been outlined in the introduction. The first “load- 
proportional” usage M.^ is based on the average load 
of country n, 



(L n ) T 

E m (L m ) 


(15) 


( 2 ) 

whereas the second “link-based” assignment Mb, asso¬ 
ciates the capacity usage ICf of a link l = (to —l n) to its 
incident nodes n and to only, 


= y. ? ■ <“) 

l(n) 


Here l(n) denotes all links l incident to node n. Both 
assignments Mn^ and Mn^ add up to ^2 n M$ = K T . 

The new flow tracing based nodal usage of the network 
can be expressed as the sum of nodal link capacities (12) 
over all links, 


M 3) = E r 


(17) 


/o\ 

Note that this assignment also adds up to Yh n ' = 
IC T . For the simplified European transmission network, 
Ain ' 1 is illustrated in Figure |6j and compared to the two 
other propositions Mn^ and Ail? from (15) and (16). 


For the three biggest countries Germany, France and 
Great Britain the new assignment (Tt| turns out to be 


smaller than the load-proportional assignment (15). For 


the next two countries Italy and Spain it is the oppo¬ 
site. This excess usage of the network becomes clearer 
when comparing Mb, directly to Mb. ■ Since France is 
central, it has large link capacities to the other four coun¬ 
tries, which explains the rather big green bar in Figure[6] 
Great Britain, Spain and also Italy are peripheral. They 
not only use the direct links to France, but also the links 
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FIG. 4: L x N matrix of the relative nodal link capacities Kj n /1C[ resulting from a 50%/50% combination of the export / 
import pictures of flow tracing. The L = 50 links are labeled on the left and the IV = 30 nodes are labeled on the bottom. 
Each row sums up to unity. 


(3) 

beyond. This explains the reduction from Mb to Mb 
for France at the expense of an increase for Great Britain, 
Spain and Italy. 

The Netherlands, Switzerland and Austria are central 

transit countries sandwiched between big countries. Con- 

( 2 ) 

sequently, their link-based assignment Mb is exception- 
ally large. The usage assignment Mb reduces these 
values and almost coincides with the load-proportional 
assignment Mn ' 1 . For Hungary, Slovakia, Czech Repub¬ 
lic, Croatia and Slovenia the situation is similar, as those 
serve as transits between numerous small Balkan coun¬ 
tries on the one hand and big Germany, Italy and France 
on the other hand. This also explains the increased blue 
bars for the Balkan countries. Peripheral countries like 
Finland, Portugal and Ireland are in a similar situation, 


leading to an excess usage of the network. 

These findings highlight that the nodal link-capacity 
( 3 ) 

assignment Mb.' based on the flow tracing network us¬ 
age measure removes the unfair allocations resulting from 
( 2 ) 

the assignment Mb. based on directly attached links. 

Compared to the load-proportional assignment Mn \ it 

(3) 

appears that to some degree M n reallocates capacity 
assignments from the big countries to medium and small¬ 
sized peripheral countries. 

V. CONCLUSION AND OUTLOOK 

Flow tracing has been developed for transmission con¬ 
gestion management and pricing in conventional power 
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FIG. 5: Relative link capacity K.f n /K.f in the export picture for the nodes DE=Germany (left) and AT=Austria (right), shown 
as color-coded links. Also shown is the average export transfer function (131, shown as color-coded nodes. 



FIG. 6: Comparison of the three capacity assignments M„ 
export and import flow tracing. 


from (151, (161 and (171, based on a 50%/50% combination of 


grids [25]. In its simplest form, it envokes the princi¬ 
ple of proportional sharing [19] , which represents a color 
diffusion on a directed flow graph embedded in the net¬ 
work. For the first time, we have applied flow tracing 
to complex renewable energy networks. Those consist of 
nodes with fluctuating renewable power generation lead¬ 
ing to fluctuating power flows on the links. The link 
capacities are determined by extreme flows and turn out 
to be quite large 03 EH- This poses the question of 
who has to pay for the link capacities. Since the flows 
are caused by the interaction between multiple injecting 
sources and ejecting sinks, the transmission capacity of 
a selected link is not only used by the two nodes directly 


attached to the link, but by more or less all nodes pro¬ 
ducing power exports and imports within the network. 
This topic has been and still is the subject of various de¬ 
bates and regulations in the European Union. For new 
infrastructure projects of common interest, the Agency 
for the Cooperation of Energy Regulators (ACER) now 
envisages the a priori cross-border allocation of the cost 
of such projects based on a cost-benefit analysis using 
detailed grid models 123 EH- With respect to exist¬ 
ing grid infrastructure, transmission system operators 
are compensated ex-post for costs induced by transient 
cross-border flows. The current mechanism, which only 
roughly estimates the influence of imports and exports 















on external grids [25] . is discussed to be replaced by a 
more appropriate one in the near future [29 j . Flow trac¬ 
ing allows to derive an expression for a fair fluctuation- 
averaged nodal usage assignment of the link capacities. 
This is showcased for a simplified European transmission 
network with a high penetration of variable renewable 
power generation, where countries are treated as nodes 
linked by interconnectors and where data from a Renew¬ 
able Energy Atlas [2j is used to determine the fluctuating 
power flows. 

This paper is intended as a first exploration of flow 
tracing applied to renewable energy networks. It will be 
interesting to extend the discussions to arbitrary renew¬ 
able penetrations ^ 1. Penetrations 0 < < 1 de¬ 

scribe intermediate stages of the transition from an elec¬ 
tricity system dominated by conventional power sources 
to one dominated by renewable resources, and flow trac¬ 
ing will allow to quantify the cost allocation of the fur¬ 
ther strengthening of the transmission grid. Furthermore 
it will be interesting to differentiate the flows between 
exporters and importers into contributions coming from 
wind power generation, solar power generation and con¬ 
ventional backup power generation. Here, flow tracing 
will help to allocate the future transmission grid exten¬ 
sions to the different forms of renewable power generation 
and to assign C0 2 emission charges directly to the con¬ 


sumers of conventional power generation in a fair manner. 

The new flow tracing based usage assignment as pre¬ 
sented in this paper acts as an ex post measure on a 
given ensemble of flow patterns in a complex renewable 
electricity network. Nevertheless, the method can also 
be used as a feedback mechanism providing incentives 
for future investments both in globally optimized scenar¬ 
ios 130] . and for individual investments and operational 
decisions in electricity markets. 

From the physics perspective, other immediate top¬ 
ics of interest are to relate flow tracing in larger net¬ 
works with significantly more nodes and links to central¬ 
ity measures of the network structure, and to study the 
transition from a global to a local coupling between ex¬ 
ports and imports once transmission capacities are con¬ 
strained. Also a comparison to other flow-based mea¬ 
sures like marginal participation, which relies on the use 
of the Power Transfer Distribution Factor (PTDF) ma¬ 
trix, will be interesting, particularly with respect to the 
role of partial counter-flows [22], [31]. Finally the relation¬ 
ship between power flows, diffusion processes and random 
walks on networks suggests to study flow tracing related 
concepts for these types of dynamics in abstract models, 
providing new insights into the fundamental properties 
of transport processes in complex systems. 
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